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Analytic  Models  for  Radiation  Induced  Loss  in 
Optical  Fibers  III.  Transient  Radiation  Effects 


1.  INTRODUCTION 

1  2 

Two  previous  HADC'  technical  reports  ’  presented  analytic  models  for  the 
losses  induced  in  optical  fibers  as  a  function  of  dose  under  steady-state  irradia¬ 
tion  conditions.  In  those  reports  it  was  assumed  that  a  state  of  at  least  quasi¬ 
equilibrium  existed  between  the  electron-hole  pairs  generated  during  irradiation 
and  the  trapping  centers  in  the  fiber.  Any  irradiation  during  which  such  equilib¬ 
rium  conditions  are  not  established  could  be  classified  as  a  transient  irradiation. 
However,  in  this  report  we  define  a  transient  irradiation  as  exposure  to  a  pulse 
of  radiation  of  short  duration  and  high  dose-rate.  The  loss  induced  in  the  fiber 
by  the  radiation  pulse  generally  decreases  (recovers)  during  the  time  following 
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the  pulse.  A  model  for  the  recovery  of  transient  induced  loss  in  optical  fibers, 
developed  to  describe  the  results  of  irradiations  performed  at  RADC*,  will  be 
presented. 

The  recovery  of  transient  radiation  induced  loss  in  optical  fibers  is  a  result 
of  the  time -dependent  recombination  of  electron-hole  pairs  generated  in  the 
fibers  by  the  radiation.  Such  recombination  processes  are  important  in  many 
phenomena  addressed  by  electronic  solid-state  physics,  such  as  the  decay  of 
luminescence  and  the  time  dependence  of  photoconductivity.  There  have  there¬ 
fore  been  a  number  of  theoretical  models  developed  to  explain  the  variety  of 
possible  recombination  processes.  Attempts  have  been  made  to  apply  these 

models  to  the  recovery  of  transient  induced  loss  in  optical  libers  (see,  for 
5  fi 

example,  Mattern  et  al  and  looney  et  al  )  with  limited  success.  It  is  not  clear, 
however,  whether  or  not  such  models  are  suitable  for  the  description  of  the  re¬ 
covery  of  induced  loss  in  optical  fibers  and  it  is  possible  that  the  phenomena 
assumed  in  developing  a  particular  model  may  not  be  those  that  occur  in  a  fiber 
even  though  the  model  may  give  a  good  fit  to  experimental  data.  Bec  ause  of 
these  uncertainties,  an  analytic  model  for  the  recovery  of  transient  radiation 
induced  loss  in  optical  fibers  was  sought  that  would  fit  the  data  well  but  not 
necessarily  be  related  to  a  physical  model. 

2.  MODKI,  SM.KCTION 

One  of  the  first  analytic  models  tested  to  describe  the  recovery  of  transient 
induced  loss  in  optical  fibers  was  a  sum  of  exponential  terms.  This  could  be 
related  to  a  possible  physical  model  in  which  each  term  represented  a  type  of 

♦In  the  transient  radiation  tests  performed  at  RADC,  20  nsec  x-ray  pulses  with 
nominal  peak  energy  of  2  MoV  and  dose  rates  to  greater  than  10^  rad/sec  were 
used.  Details  of  the  irradiations  are  described  elsewhere.  3.  4 

3.  Wall,  J.A.,  Rosen,  H.  ,  and  Jaeger,  H.  (1981)  Temperature  response  of 

germanium  phosphosilicate  optical  fibers  under  irradiation,  in  Physics  of 
Fiber  Optics,  Advances  in  Ceramics,  U.  Bcndow  and  S.  S.  Mitre.,  Fds.  , 
American  Ceramic  Society  2:398. 
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trapping  center  with  a  unique  temporal  probability  for  the  release  and  subsequent 
recombination  of  an  electron  or  hole.  It  was  found  that  most  of  the  RADC  data  on 
transient  induced  loss  in  optical  fibers  could  be  fit  to  this  model  using  three  ex¬ 
ponential  terms.  However,  another  model  was  sought,  because  the  exponential 
model  was  difficult  to  fit  to  the  data  since  it  required  the  determination  of  six 
constants  (two  for  each  term).  Also,  its  usefulness  in  providing  a  physical  inter¬ 
pretation  of  the  data  was  questionable  because  any  continuous  mathematical  func¬ 
tion  can  be  represented  by  the  sum  of  a  sufficient  number  of  exponential  terms. 

Since  a  number  of  the  theoretical  models  for  electron-hole  recombination 
processes  could  be  approximated  by  a  function  of  the  form  cfn,  where  t  is  the 
time  and  c  and  n  are  constants,  a  model  of  this  form  was  tested.  It  was  found  to 
fit  the  data  fairly  well  over  much  of  the  time  span  covered  by  the  data,  but  has 
the  obvious  disadvantage  of  tending  to  infinity  at  short  times.  To  avoid  this  dis¬ 
advantage  the  function  was  arbitrarily  modified  to  the  form: 


L  =  radiation  induced  loss 
t  =  time  after  radiation  pulse,  and 
a,  b,  and  n  are  constants. 

This  expression  was  found  to  fit  the  transient  data  quite  well.  Although  not  re¬ 
lated  to  a  physical  model,  Eq.  (1)  has  the  form  theoretically  derived  for  the  bi- 
molecular  recombination  process  if  n  *  1.  Also,  according  to  Haug,  it  could 
be  an  approximation  for  thermally  activated  recombination  and  according  to 

O 

Drown,  it  could  represent  recombination  from  an  exponential  distribution  of 
traps. 

To  give  a  complete  description  of  transient  radiation  induced  loss  data  a 
constant  or  term  that  decreases  slowly  with  time  should  be  added  to  Eq.  (1)  to 
account  for  permanent  or  long-term  induced  loss.  However,  because  the  total 
dose  received  by  the  optical  fibers  during  the  transient  radiation  tests  at  RADC 
was  ho  small  (lews  than  :•()  rnd),  this  factor  will  be  neglected. 


7.  lluug,  A.  (11172)  Theotcticnl  Solid  State  Physics  Volume  2,  Pergamon  Press, 

pp. 

ti.  Ilrowii,  K.C.  (ll'M.)  The  Physics  of  Solids,  W.A.  It.  njamin  Inc.  ,  N.Y., 
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3.  DATA  FITTING  PROCEDURE 


To  determine  the  constants  in  Eq.  (1)  for  a  specific-set  of  transient  induced 
loss  data,  it  is  convenient  to  express  the  data  in  terms  of  reciprocal  loss  as  a 
function  of  time.  Instead  of  Eq.  (1),  we  then  fit 


L 


-1 


(2) 


and  determine  the  reciprocal  of  a,  b/a,  and  n.  This  expression  can  be  readily  fit 
using  the  "Gauss -Newton"  method  mentioned  in  Ref.  2.  However,  good  fits  have 
been  obtained  using  the  following  procedure,  which  is  also  useful  for  obtaining 
estimates  of  the  constants  needed  for  the  Gauss -Newton  method. 

The  first  derivative  of  Eq.  (2)  is 


_d 

dt 


a/1) 


i 


(3) 


Using  a  finite  difference  method  to  approximate  the  differential  from  the  data  and 
a  standard  procedure  for  fitting  a  power  function  (see  Ref.  1  for  example)  a  value 
for  n  (actually  n-1)  is  obtained.  When  the  time  data  are  converted  to  the  nth  pow  er 
of  time,  a  linear  fit  to  the  reciprocal  of  the  induced  loss  [  Kq.  (2))  gives  the 
required  values  of  l/a  and  b/a. 

It  should  be  noted  that  when  the  data  are  fitted  to  the  approximation  of  Eq.  (3), 
low  values  of  the  correlation  coefficient  for  the  fitting  procedure  may  frequently 
be  obtained.  This  is  probably  due  to  the  approximation  of  the  differential  and 
seems  to  have  little  relation  to  the  final  result,  since  high  values  of  the  correla¬ 
tion  coefficient  have  usually  been  obtained  for  the  fit  to  Eq.  (2). 
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4.  RESULTS 


tr.y 


Figure  1  shows  the  results  of  a  fit  of  Eq.  (1)  to  data  for  the  recovery  of 
transient  radiation  induced  loss  in  a  germanium  -boros ilieate  core  optical  fiber. 
The  fit  is  very  good  and  the  deviation  of  the  calculated  curve  from  the  data  points 
is  generally  well  within  the  estimated  ±  3  dli/km  precision  of  the  data.  Figure  2 
shows  two  fits  of  Eq.  (1)  to  transient  data  for  a  germanium -phosphosilicate  core 
optical  fiber.  The  first  fit  (solid  curve)  utilized  all  of  the  data  up  to  100  usee.  As 
can  be  seen  in  the  figure  the  fit  is  not  very  good  below  2  usee.  For  the  second  fit 
(dashed  curve),  points  taken  near  the  ends  of  the  oscilloscope  traces  (the  data 
were  read  from  four  separate  oscilloscope"  beams),  whi<  h  are  difficult  to  read 
with  reasonable  accuracy,  were  eliminated.  The  second  fit  is  clearly  much  better 
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Figure  1.  Fit  of  Eq.  (1)  to  Data  for  the  Recovery  of  Transient  Radiation  In¬ 
duced  Loss  in  a  Germanium -boros ilicate  Core  Optical  Fiber.  The  irradia¬ 
tion  conditions  were:  pulse  width,  20  nsec;  dose-rate,  1.  3  X  109  rad/sec; 
peak  x-ray  energy,  2  MeV;  temperature  25C;  source  wavelength  820  nm 

than  the  first.  In  both  figures,  the  fitted  curves  deviate  somewhat  more  fri.m  the 
data  for  times  longer  than  100  p sec.  This  could  be  due  to  small  differences  in 
the  calibrations  of  the  oscilloscope  amplifiers  and/or  time  bases,  which  may  in¬ 
fluence  the  accuracy  of  the  data  in  this  time  range  where  the  signal  to  noise  ratio 
is  low.  This  possibility  is  most  evident  in  Figure  2  where  the  data  shows  a 
greater  loss  remaining  at  400  psec  than  at  100  psec,  a  very  unlikely  event. 

Transient  data  for  more  than  12  optical  fibers  were  fit  to  Eq.  (1)  with  re¬ 
sults  similar  to  those  shown  in  Figures  1  and  2.  Good  fits  were  also  obtained  on 
fibers  irradiated  at  temperatures  of  -55  and  +125C.  Considering  the  time  spans 
involved  as  well  as  the  inherent  "noise"  associated  with  the  acquisition  of  tran¬ 
sient  data,  this  is  remarkable  for  an  expression  requiring  the  determination  of 
only  three  constants. 


5.  DISCUSSION 

The  success  of  Eq.  (1)  in  describing  the  recovery  of  transient  radiation  in 
duci  d  loss  in  optical  fibers  runimt  be  explained  except  that  il  may  be  n  good 
approximation  to  several  physical  models  developed  to  explain  electron-hole 
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Figure  2.  Fit  of  Eq.  (1)  to  Data  for  the  Recovery  of  Transient  Radiation  In¬ 
duced  Loss  in  a  Germanium-Phosphosilicate  Core  Optical  Fiber.  Fit  1  (solid 
curve),  all  data  used.  Fit  2  (dashed  curve),  data  of  questionable  accuracy 
eliminated  (see  text).  The  irradiation  conditions  were:  pulse  width,  20  nsec; 
dose  rate,  1.  1  X  10^  rad/sec;  peak  x-ray  energy,  2  MeV;  temperature,  25C; 
source  wavelength,  820  nm 

recombination  processes  in  other  materials.  Evidence  that  Eq.  (1)  has  some  de¬ 
gree  of  "universality"  is  given  by  the  fact  that,  with  the  addition  of  a  constant  to 
account  for  permanent  induced  loss,  it  could  be  used  to  fit  data  on  the  recovery 
of  induced  loss  in  optical  fillers  following  steady-state  irradiation.  [Specifically, 
Eq.  (1)  plus  a  constant  gave  very  good  fits  to  the  data  shown  in  Figures  3  and  4 
of  Ref.  4.  ]  However,  there  is  insufficient  information  available  at  this  time  to 
determine  whether  or  not  the  constants  obtained  in  fitting  Eq.  (1)  to  data  may  be 
dose -rate  and/or  dose  dependent.  Any  extrapolation  of  the  results  of  data  fits  to 
Eq.  (1)  should  therefore  be  performed  with  this  uncertainty  in  mind.  Also, 
because  of  the  complexity  of  the  dynamics  of  electron-hole  trapping  and  release 
processes,  it  is  possible  that  Eq.  (1)  may  not  he  applicable  to  data  obtained  for 
irradiations  performed  at  much  higher  doses  and  dose-rates  or  taken  at  shorter 
times  than  those  for  which  the  data  fits  described  in  the  report  were  performed. 
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§  MISSION  I 

^  Rome  Air  Development  Center  % 

n  RA VC  plam  and  executes  research ,  de.vzlopme.nt,  test  v 

•  and  selected  acquis Ition  programs  in  support  of  V 

&  Command,  Control,  Communications  and  Intelligence 

*t  (C3 1)  activities.  Technical  and  engineering  £ 

&  support  within  areas  competence  is  provided  to  § 
D  CSV  Program  Offices  { POs )  and  other  ESP  elements  x 

?  to  perform  elective  acquisition  o  f  C$  I  systems.  § 

q  The  areas  of  technical  competence  inclucfe 
\  communications,  command  and  control,  battle  Q 

&  management,  information  processing ,  surveillance 
sensors,  intelligence  data  collection  and  handling, 

£  solid  state  sciences ,  electromagnetics ,  and 
n  propagation,  and  electronic,  maintainability, 

*  and  compatibility. 
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